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Abstract 



Hypcrfine induced ls2s 1 So — > Is 2 1 Sq Ml transition probabilities of He-like ions have been calcu- 
lated from rclativistic configuration interaction wavefunctions including the frequency independent 
Brcit interaction and QED effects. Present results for 151 Eu and 155 Gd are in good agreement with 
previous calculations [Phys. Rev. A 63, 054105 (2001)]. Electronic data are given in terms of a 
general scaling law in Z that, given isotopic nuclear spin and magnetic moment, allows hypcrfine 
induced decay rates to be estimated for any isotope. The results should be helpful for future 
experimental investigations on QED and parity non-conservation effects. 



PACS: 31.30.Gs, 31.15.V-, 31.15.ag 



The hypcrfine structure of atomic levels is caused by the interaction between the electrons 
and the non-central electromagnetic multipoles of the nucleus. The interaction, although weak, 
not only shifts and splits the individual J levels, but also mixes wavefunctions with different J 
quantum numbers pQ . The mixing may open new decay channels and induce radiative transitions 
such as J = — > J' = 0. The study of hyperfine induced transitions (HIT) is attractive in view of 
getting a more subtle understanding of electron correlation and relativistic effects [2j [3] [4] , testing 
the standard model [5], obtaining nuclear properties [6l [7j |8], developing ultra-precise atomic 
clocks [TOl HU Q~2j , analyzing stellar spectra [13] , and determining isotopic abundance ratios and 
electron densities of plasma [33] [T3] . 

Hc-like systems are interesting due to their relative simplicity, and important effects such as 
electron correlation [15], quantum electrodynamics (QED) [TTlfTSlfTO] and parity non-conservation 
(PNC) phenomena [30] dTJ HH H3] can be investigated in detail. Gorshkov and Labzovskii 20] 
and Labzowsky et al. [23| have proposed that the mixed hyperfine- and weak-quenching can be 
used to test parity- violation effects. The one photon transition ls2s 1 S'o — > Is 2 1 5o of Hc-like ions 
is considered a good candidate for these tests and experiments will be carried out at GSI [24] . 
Therefore, accurate hyperfine induced ls2s 1 S'o — > Is 2 1 5'o Ml transition probabilities of Hc-like 
ions are important. Relevant data are however still insufficient and in response to this we have 
performed systematic calculations along the He-like iso-electronic sequence using GRASP2K [25] 
based on the multi-configuration Dirac-Hartree-Fock method and the HFST [25] package. From 
computed hyperfine induced transitions for 14 ions and corresponding electronic quantities we 
have derived a scaling formula in Z that allows the induced transition rate to be computed for 
any isotope. 

When the hyperfine interaction is included the wavefunctions of the combined electronic and 
nuclear system can be represented as 

\YFM F ) = Y J c i \liUFM F ). (1) 

i 

The zero-order functions \jIJFMp) in the expansion are coupled products of electronic \~yJMj) 
and nuclear \IMj) wavefunctions. The Is 2 1 Sq ground state is well represented by a single term. 
For ls2s 1 Sq only the interaction with ls2s 3 Si is important (see Figure 1) and the wavefunction 
can be approximated by the expansion 

"ls2s ^oIF") =c \ls2s 1 S Q IF) + C\\ls2s 3 S 1 IF), (2) 

where / is the nuclear spin and F(= I) the total angular momentum quantum number. Magnetic 
quantum numbers are suppressed for brevity. The use of quotation marks in the left-hand wave- 
function emphasizes the fact that the notation is just a label indicating the dominant character 
of the eigenvector. The mixing coefficient c\ is obtained in first order perturbation theory as the 
ratio between the hyperfine matrix element and the unperturbed energy differences 

(lfga 3 Si IF\H hfs \ls2s 1 S IF) 
Cl E(iS ) - E^S,) ■ [d) 

The one-photon ls2s 1 5o — » Is 2 1 S'o Ml transition becomes allowed via the mixing of ls2s 3 S\ 
and the decay rate in s _1 is given by 

A nolo 12 1cm 2.69735 x 10 13 c 2 5mi 

A HIT (ls2s S -> Is So) = ^3 , (4) 

where 

Smi = (Is 2 1 5 ||M«||1 S 2 S 3 5 1 ) ' (5) 



is the line strength computed from unperturbed wavefunctions. A is the wavelength in A for the 
transition. The reader is referred to Q3] HZ] [55] for details of the derivation. 



The electronic wavefunctions were computed using the GRASP2K program package [35] ■ Here 
the wavefunction for a state labeled 7 J is approximated by an expansion over jj-coupled config- 
uration state functions (CSFs) 

\ 7 J} = J2di$i. (6) 

z 

In the multi-configuration self-consistent field (SCF) procedure both the radial parts of the or- 
bitals and the expansion coefficients are optimized to self-consistency. In the present work a 
Dirac-Coulomb Hamiltonian was used with the nucleus described by an extended Fermi charge 
distribution [5^. The multi-configuration SCF calculations were followed by relativistic CI calcu- 
lations including the frequency independent Breit interaction and leading QED effects. 

For the low charged ions the main uncertainties in the calculation come from electron corre- 
lation effects. To build a reasonable correlation model and control the accuracy we performed 
tentative calculations of transition energies and the (ls2s 3 Si IF\Hhf s \ls2s 1 SoIF) off-diagonal 
hyperfme matrix element in 12 C. In the calculations the wavefunctions for Is 2 So, ls2s 
ls2s 3 Si were determined simultaneously in extended optimal level (EOL) calculations [30]. All 
CSFs that could be built from an active set of orbitals were included in the expansion. The 
orbital set was systematically increased by adding layers of new orbitals. The largest active set 
included relativistic orbitals with principal quantum number n ^ 7. Due to stability problems in 
the relativistic SCF procedure only the outermost layers of orbitals could be optimized each time. 
The frequency independent Breit interaction and leading QED effects were added in subsequent 
relativistic CI calculations. The results for 12 C are shown in Table 1. The first column in this 
table represent the active set of orbitals involved in each step of the calculation. As can be seen 
from the table, five energy-optimized layers of orbitals are needed to converge the off-diagonal 
hyperfme interaction matrix element between ls2s 3 Si and ls2s ^^Sq at a sub per mill level. 

Based on the above analysis, we performed calculations for other He-like ions using an active 
orbital set with n ^ 5. The Breit interaction and main QED corrections were included. Values of 
nuclear magnetic dipole moments for the different isotopes were adopted from the compilation by 
Stone |32j . In Table 2 we display the off-diagonal hyperfme matrix elements and corresponding 
mixing coefficients c\. The given values are not corrected for the the anomalous magnetic moment. 
For the mixing coefficient of 151 Eu the magnitude is agreement with the value of [53], but with 
a difference in sign due to different definitions of the phase factor in the hyperfme interaction 
matrix element. The difference in sign does not influence the final hyperfme induced transition 
probability. Wavelengths A and line strengths 5m 1 for the ls2s 3 Si — Is 2 1 So Ml transition needed 
for the quenching rate are taken from accurate relativistic CI calculations by Johnson et al. [18] . 
The values are presented in Table 3 for the convenience of the reader. The hyperfme induced 
Is2s 1 So — > Is 2 1 So Ml transition rate and corresponding wavelengths of He-like ions are given 
Table 4. Previous theoretical results of wavelength [33] are compared with present calculations 
in this table. The agreement between our value for 151 Eu and previous theoretical values by 
Labzowsky et al. |23j is very good. 

To predict the transition rate for any isotope in the iso-electronic sequence we follow Brage et 
al. [14] and factorize the hyperfme induced transition rate into nuclear and electronic parts 

A HIT (ls2s l S -> Is 2 l S ) = + l/I)A el (ls2s l S -> Is 2 l S ). (7) 

The electronic part A el has a smooth behavior along the iso-electronic sequence making interpo- 
lation possible. From the data in Table 4 we obtain a fit of the form 

A el = 1.9728 x io- 19 Z 14 065 (8) 

where Z is the atomic number. The fit is shown in Figure 2. In order to show clearly the trend of 
A el with Z, log(A) was plotted in this picture. Using the fitting formula we estimate the probability 
of 155 Gd with nuclear spin / = 3/2 and nuclear dipole moment /j = — 0.2591/xjv to 5.60 xlO 5 s _1 . 
This value is in good agreement with the theoretical value 5.8 xlO 5 s _1 given in [53]. The fitting 
formula (8) is expected to provide accurate values, but further examination through experiment 
and theory is still needed. 



To sum up we have calculated hypcrfinc induced ls2s 1 Sq — > Is 2 1 Sq Ml transition probabilities 
of He-like ions using GRASP2K j^S] based on multi-configuration Dirac-Fock method and the 
HFST [HI] package. Electron correlation effects were included in a systematic way. The Brcit 
interaction and QED effects were included in subsequent relativistic CI calculations. A scaling 
law in Z was derived for the electronic quantities. The scaling law allows hyperfine induced 
transition probabilities to be estimated for any isotope. 
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Figure 1: Level structure of the Is 2 , ls2s and ls3s configurations in Hc-likc C and Ir. 

Table 1: Transition energies in cm -1 and off-diagonal hyperfine interaction matrix element in 
a.u. for 12 C from calculations with increasing active sets. AEi = E(ls2s 3 Si) — E(ls 2 1 5'o), 



AE 2 = E(ls2s l S ) 


-E(ls2s 


A£ 3 


= E(ls2 


s l S )-E(ls2s 3 


active set 


AE 1 


AE 2 


AE 3 


Matrix Element 


2slp 


2407416 


2451264 


43847 


4.3613[-6] 


3s2pld 


2410850 


2455165 


44315 


4.5773[-6] 


4s3p2dlf 


2411282 


2455181 


43899 


4.5641 [-6] 


5s4p3d2flg 


2411308 


2455148 


43840 


4.5576 [-6] 


6s5p4d3f2glh 


2411422 


2455224 


43802 


4.5576 [-6] 


7s6p5d4f3g2hli 


2411475 


2455265 


43791 


4.5575[-6] 


QED corrections 


2411439 


2455237 


43798 


4.5572[-6] 


Experiment 1311 


2411271 


2455026 


43755 





Table 2: Off-diagonal hyperfine interaction matrix elements (a.u.) and hyperfine mixing coeffi- 
cients ci of Hc-likc ions. Nuclear magnetic dipolc moment fij arc from |32j . 



Mixing Coefficient 



Isotope 


Z 


I 


M/ 


Matrix Element 


This work 


Rcf. 23^ 


ia c 


6 


1/2 


0.7042 


4.5576[-6] 


2.2817[-5] 




19 F 


9 


1/2 


2.628868 


5.7260[-5] 


1.6993[-4] 




28 Si 


14 


1/2 


-0.55529 


-4.5667[-5] 


-8.1505[-5] 




47 Ti 


22 


5/2 


-0.78848 


-1.7489 [-4] 


-1.8463[-4] 




57 Fe 


26 


1/2 


2.56277 


4.8981[-5] 


4.2385 [-5] 




70 Ga 


31 


3/2 


2.6289 


1.7944 [-3] 


1.2473[-3] 




85 Rb 


37 


5/2 


1.35298 


1.5213[-3] 


8.4420[-4] 




97 Mo 


42 


5/2 


-0.9335 


-1.5826[-3] 


-7.3716[-4] 




103 Rh 


45 


1/2 


-0.884 


-2. 7542 [-3] 


-1.1627[-3] 




117 Sn 


50 


1/2 


-1.00104 


-4.4439 [-3] 


-1.6030[-3] 




131 Xc 


54 


3/2 


0.6908 


2.9782[-3] 


9.5300[-4] 




151 Eu 


63 


5/2 


3.4717 


2.3846[-2] 


5.8619[-3] 


-5.83[-3] 


175 Lu 


71 


7/2 


2.2323 


2.3187[-2] 


4.5372 [-3] 




193 j r 


77 


3/2 


0.1637 


2.6895[-3] 


4.4400[-4] 





Tabic 3: Wavelengths A in A, rates Ami in s 1 and line strengths Smi in a.u. for the ls2s 3 Si — 
Is 2 1 Sq Ml transition of He-like ions. From Johnson et al. [18] . 



Ion 


Z 


A 


Ami 


Smi 




C 


6 


41.470 


4.860[1] 


3.8550| 


-7} 


F 


9 


17.152 


3.621[3] 


2.0323| 


-6] 


Si 


14 


6.7402 


3.598[5] 


1.2253| 


-5] 


Ti 


22 


2.6368 


3.750[7] 


7.6467| 


-5] 


Fe 


2G 


1.8682 


2.075[8] 


1.5046| 


-4] 


Ga 


31 


1.3002 


1.257[9] 


3.0719| 


-4] 


Rb 


37 


0.90272 


7.714[9] 


6.3110| 


-4] 


Mo 


42 


0.69459 


2.842[10] 


1.0592| 


-3] 


Rh 


45 


0.60198 


5.792[10] 


1.4054| 


-3] 


Sn 


50 


0.48342 


1.726[11] 


2.1681 


-3] 


Xe 


54 


0.41151 


3.846[11] 


2.9807| 


-3] 


Eu 


63 


0.29715 


1.943[12] 


5.6695| 


-3] 


Lu 


71 


0.22993 


6.950[12] 


9.3969| 


-3] 


Ir 


77 


0.19266 


1.672[13] 


1.3300| 


-2] 



Table 4: Hypcrfine induced rates Ah it in s 1 and corresponding wavelength A in A for the 

ls2s 1 5o — > Is 2 1 Sq transition in He-like ions. 

A(A) A H it 

Isotope Z I m This work Rcf. 33J This work Rcf.[23j 

'^G 6 T72 0.7042 40.7306 40.7304 2.6704[-8] 

19 F 9 1/2 2.628868 16.9361 16.9404 1.0862 [-4] 

28 Si 14 1/2 -0.55529 6.6835 6.6848 2.4515[-3] 

47 Ti 22 5/2 -0.78848 2.6214 2.6225 1.3010[0] 

57 Fe 26 1/2 2.56277 1.8584 1.8594 3.7865[-l] 

70 Ga 31 3/2 2.6289 1.2940 1.9830[3] 

85 Rb 37 5/2 1.35298 0.8987 5.5705[3] 

97 Mo 42 5/2 -0.9335 0.6916 0.6923 1.5643[4] 

103 Rh 45 1/2 -0.884 0.5994 7.9312[4] 

117 Sn 50 1/2 -1.00104 0.4814 0.4820 4.4904[5] 

131 Xe 54 3/2 0.6908 0.4098 0.4104 3.5374[5] 

151 Eu 63 5/2 3.4717 0.2958 6.7643[7] 6.8[7] 

175 Lu 71 7/2 2.2323 0.2289 1.4508[8] 

193 Ir 77 3/2 0.1637 0.1917 3.3463[6] 



